Conductance of single-atom magnetic junctions: A first-principles study 
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We present a first-principles investigation to show that the contact conductance of a half conductance quantum 
(Go/2) found previously does not generally hold for single-atom magnetic junctions composed of a tip and 
an adatom adsorbed on a surface. The contact conductance of the Ni-Co/Co(l 1 1) junction is approximately 
00/2, while for the Co-Co/Co(l 1 1), Ni-Ni/Ni(l 1 1), and Ni-Ni/Ni(001) junctions the contact conductances are 
O.SOGo, I.55G0, and 1.77Go, respectively. The deviation from Go/2 is mainly caused by the variation of 
the spin-down conductance largely determined by the minority d orbitals, as the spin-up one changes little for 
different junctions. 
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The electrical transport property of the atomic contact has 
attracted much attention for its enormous potential in future 
atomicscale electronic devices. The conductance of a spin- 
degenerated atomic-size contact is quantized in units of Go = 
2e^/h, where e is the proton charge and h is the Planck's con- 
stant according to the Landauer's formula|l]. In the case of 
magnetic systems the spin degeneracy is removed, and each 
spin-polarized channel can contribute up to Go/ 2 to the total 
conductance. 

The conductance of atomic junctions has been extensively 
studied by experimental techniques such as the scanning 
tunneling microscope (STM) [2-^, mechanically controlled 
break junctio ns ( MBCJ) technique [7-11], and electrochem- 
ical methods ifl^, ITsll . Conductances of fractional and inte- 
ger Go have been reported for both nonmagnetic and mag- 
netic junctions. Of the particular interest is the observation of 
Go/2 conductance, e.g., for the noble Pt, Pd, and magnetic 
Co sing le-atom chains f§\, as well as Cu and Ni nanowires 
lfT2lll4ll . although the controversy still exists by arguing that 
the contamination of H2 may cause the Go/2 conductance 
ifToll . Theoretically, the transport property of atomic junctions 
has been investigated intensively by using either tight-binding 
methods fT5'-'l7] or density functional theory (DFT) com- 
bined with quantum transport calculations based on nonequi- 
librium Green's function (NEGF) formalism [18-27]. How- 
ever, theoretical calculations have hardly obtained the Go/ 2 
conductance so far. This may be attributed to the difference 
between the simulation model and the real junction structure 
which is usually difficult to know experimentally. 

Recently, a STM experiment found that in the contact 
regime the conductance is Go/2 for a single-atom magnetic 
junction composed of a feiTomagnetic Ni tip in contact with 
a single Co adatom on a ferromagnetic Co island jljtj. This 
experiment provides a more detailed atomic geometry of the 
junction than those ever reported, which enables a precise 
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theoretical modeling to explore the origin of the Go/2 con- 
ductance. Calculations show that the conductance of Go/2 
is not due to a fully polarized single channel but a combi- 
nation of the partially open majority channels and the sup- 
pressed minority channels [24]. However, it still remains un- 
known whether the Go/2 conductance exists generally for 
other single-atom magnetic junctions, and how factors like 
junction structure and species affect the conclusion. In this 
work, we perform first-principles calculations to investigate 
this issue and find that the Go/2 conductance is not a gen- 
eral behavior for several other single-atom magnetic junc- 
tions. This is because the minority-spin conductance is sensi- 
tive to the junction species. 

We first study the magnetic junction fabricated by Neel et 
al.[3] The system is modeled by a tip-adatom/surface junction 
(denoted by Ni-Co/Co(l 1 1)) as displayed in FiglTJa). The sin- 
gle Co adatom is adsorbed on the fee site of the fee Co(l 1 1) 
surface represented by a three-layer slab, with each layer con- 
taining 4x4 atoms. The Ni tip is modeled by a single Ni apex 
atom and a Ni(lll) monolayer on a two-layer Co(lll) slab 
(test calculations using three Ni(lll) layers on a two-layer 
Co(l 1 1) slab show only minor effect on the result). The tip- 
apex atom is placed above the Co adatom in the z direction. 
In transport calculations, these atoms construct the scattering 
region, and three additional Co( 111) layers are added at the 
two ends of the scattering region, respectively, to mimic the 
left and right electrodes (leads). For the other three junctions 
considered in this work, i.e., Co-Co/Co(l 1 1), Ni-Ni/Ni(l 1 1), 
and Ni-Ni/Ni(001), the treatment is similar except that a four- 
layer slab is used for the Ni(OOl) surface and four additional 
atomic layers are added to represent the leads of the Ni(OOl) 
junction, as illustrated in Fig (Tib). 

In this work, the ferromagnetic spin configuration is 
adopted for the whole junction. Our calculations show that if 
the junction (including the tip and the substrate) is considered 
as one system, the ferromagnetic configuration is the ground 
state, as was also found in a previous calculation for the Ni- 
Co/Co(l 1 1) junction |22]. If the junction is considered as two 
separated system (i.e., the tip and the substrate) the antifeiTO- 
magnetic configuration between them may be also possible. 
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FIG. 1 : Model system for the tansport calculation of the single-atom 
magnetic junction, (a) The Ni-Co/Co(l 11) junction: A Ni tip (green 
balls) is placed above the Co adatom on the Co(ll 1) surface, (b) The 
Ni-Ni/Ni(001) junction. The transport is along the Z direction. 



as considered in Ref. However, the only consideration 

of the ferromagnetic configuration does not affect the pur- 
pose and conclusion of this work: The contact conductance 
of 6*0/2 is not a general behavior for different single-atom 
magnetic junctions. 

The scattering region is optimized using VASP code llZTl . 
The bottom layer of the surface and the top layer of the tip are 
fixed during the structure optimization while the other atoms 
are fully relaxed until the maximum force is smaller than 
O.OleV/A. Projector augmented-wave method [28] is used for 
the wave function expansion with an energy cutoff of 300eV. 
The PW91 version [29] of the generalized gradient approxi- 
mation (GGA) is adopted for the electron exchange and cor- 
relation. 

For the quantum transport calculation, we adopt the NEGF- 
DFT approach ll30l [31I1 which combines the NEGF formula 
for transport with ab initio DFT calculation for electronic 
structure. In practice, the infinitely large open system is di- 
vided into three parts: left lead, right lead, and scattering re- 
gion, as mentioned above. The self-consistent Kohn-Sham 
Hamiltonian of the scattering region and the self-energies of 
two semi-infinite leads are used to construct a single-particle 
Green's function from which the transmission coefficient (T) 
at any energy is calculated. The conductance G then follows 
from a Landauer-type relation. For the DFT part, we use a 
numerical basis set to expand the wavefunction [35[: A dou- 
ble zeta plus polarization basis set (DZP) is adopted for all 
atomic species. The Perdew-Burke-Ernzerhof (PBE) jssll ver- 
sion of GGA is used for the electron exchange and correlation 
and the optimized Troullier-Martins pseudopotentials [ 34] are 
used for the atomic cores. We define the spin-polarization ra- 
tio at Fermi energy as P = (Tf - T;)/(r-f- + T^), where 
and denote the transmission coefficient of the majority and 
minority spin, respectively. 

The conductance of the magnetic Ni-Co/Co( 111) junction 
is given in Fig|2|a) as a function of the tip height (i.e., the dis- 
tance between the tip-apex atom and the surface before the re- 
laxation). It shows that as the tip height decreases, the conduc- 




5.2 5.4 
Tip Height (A) 



FIG. 2: (a) The spin-dependent conductance of the Ni-Co/Co(lll) 
junction, and (b) the adsorption height of the Co adatom relative to 
the surface for different tip heights. 



tance increases and shows a faster change in the transition re- 
gion around 5. 4-5. 2 A, and then increases slowly in the contact 
region (below 5. 2 A). The conductance data in the transition 
and contact regions can be approximated by two straight lines, 
and their intersection point defines the contact conductance, 
as used in Refs.[3, 35]. According to this definition, we ob- 
tain a contact conductance of 0.48 Gq for the Ni-Co/Co(lll) 
junction, which agrees very well with the experimental result 
[3]. The two spin components of the conductance give a spin- 
polarization ratio of 0.46 for the tip height of 5. 2 A, which 
is also in good agreement with the previous first-principles 
calculation[24]. 

Figl2b) gives the adatom's adsorption height for different 
tip heights. On can see that in the contact region the adsorp- 
tion height is obviously larger than that in the tunneling re- 
gion, and so is the conductance. This is a result of the en- 
hanced adatom-tip interaction in the contact region. An im- 
portant thing to note is that the adsorption height experiences 
a large jump for the tip height around 5.3A, where the con- 
ductance also rapidly increases as shown in Fig|2|a). This 
indicates that the adsorption height of the adatom has a direct 
influence on the conductance because the resulting adatom- 
tip separation determines the coupling strength (the degree of 
orbital overlap). Physically, the adsorption height and the re- 
sulting adatom-tip separation for a specific tip height is deter- 
mined by the interaction between the tip and the substrate, 
which may depend on the junction structure and species. 
Therefore, to investigate the influence of these two factors on 
the conductance is important to find out whether a Go /2 con- 
ductance is a general result for other magnetic single-atom 
junctions. 

Next, we investigate the conductance of the following 
three junctions, Co-Co/Co(lll), Ni-Ni/Ni(lll), and Ni- 
Ni/Ni(001). The first two have a similar structure as Ni- 
Co/Co(lll) but are of different species. Thus by compar- 
ing their results, we can evaluate the role of species in de- 
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FIG. 3: (a) Conductance of the Co-Co/Co(l 1 1), Ni-Ni/Ni(l 1 1), and 
Ni-Ni/Ni(001) junctions, respectively, (b) and (c) show the spin-up 
and spin-down conductances of these three junctions and of the Ni- 
Co/Co(l 1 1) junction, respectively. 



termining the conductance of the magnetic junctions. The 
Ni-Ni/Ni(001) and Ni-Ni/Ni( 111) junctions are of the same 
species but have different structures, which enables us to ex- 
plore the influence of the atomic structure on the conductance. 

Fig 13 a) displays the conductance as a function of the tip 
height, and FigslSb) and (c) give the two spin components of 
the conductance, for the Co-Co/Co(lll), Ni-Ni/Ni(lll) and 
Ni-Ni/Ni(001) junctions, respectively. We find that their con- 
tact conductances are O.SOGo, 1.55Go, and 1. 11 Go with the 
corresponding spin-polarization ratios of 0.1, -0.45, and -0.33, 
respectively. Obviously, the contact conductances of the three 
junctions deviate largely from Gq/2 and also differ much from 
each other. 

Specifically, the spin-up conductance of the Ni- 
Co/Co( 111), Co-Co/Co( 111), and Ni-Ni/Ni( 111) junctions 
are very close (see FigHb)), while their spin-down conduc- 
tances differ from each other significantly (see Figl^c)). 
Since these three junctions are of the similar structure, this 
difference is due to the different junction species which 
give rise to different local electronic states and tip-adatom 
couplings. 

We find that the spin-up conductance is mainly contributed 
by the majority s, Pz, and dz (sum of dz^^r^, dxz and dyz 
with z being the transport direction) orbitals and their contri- 
butions to the majority conductance are almost the same for 
the four different magnetic junctions |[24ll . Consequently, the 



FIG. 4: (a) The majority and (b) minority PDOS of the Co adatom 
in the Co-Co/Co(l 1 1) junction and that of the Ni adatom in the Ni- 
Ni/Ni( 11 1) junction with a tip height of 4.4A. (c) The minority s, pz 
and dz PDOS of the Co adatom in the Co-Co/Co(ll 1) junction with 
4.4Atip height. 



spin-up conductance is little changed for the different junc- 
tions. As an example, we give in FigUa) the majority-orbital- 
projected density of states (PDOS) of the Ni adatom in the 
Ni-Ni/Ni(l 1 1) junction and that of the Co adatom in the Co- 
Co/Co(l 1 1) junction with a tip height of 4.4A. It can be seen 
that the majority PDOS at the Fermi energy are very close 
for the two adatoms, leading to the similar spin-up conduc- 
tances for the two junctions (see Fig[3jb)). In contrast, the 
spin-down conductance is basically determined by the minor- 
ity dz orbitals. Their contribution to the minority PDOS of 
the adatom at the Fermi energy is much larger than the s and 
Pz orbitals, as shown in Fig|4jc) for the Co-Co/Co(l 1 1) junc- 
tion. Importantly, as displayed in FigHfb), the minority dz 
PDOS of the adatom at the Fermi energy varies significantly 
between the junctions, leading to the very different spin-down 
conductances for these two single-atom magnetic junctions. 

Now let us have a look at the effect of the junction struc- 
ture. By comparing the results for the Ni-Ni/Ni(l 1 1) and Ni- 
Ni/Ni(001) junctions (see Fig©, we find that both the spin- 
down and spin-up conductances are different, causing a 0.2Go 
difference in the total conductance as shown in Figl3ja). 
However, this difference is smaller than those between the 
Co-Co/Co( 111), Ni-Co/Co( 1 1 1 ), and Ni-Ni/Ni( 111) junctions 
caused by the different junction species (see Figl3]l. Our re- 
sults show that the Go/2 conductance does not generally hold 
for the single-atom magnetic junctions investigated in this 
work. This is due to the the large variation of the spin-down 
conductance which largely depends on the junction species as 
well as the junction structure. 

Finally, we would like to comment on the major contribu- 
tion from the majority or the minority to the total conduc- 
tance. In Ref.0], the authors investigated the conductance 
of single-atom magnetic junctions of Cr-Co/Fe(l 10) and Cr- 
Cr/Fe(110). It was found that the conductance is mainly de- 
termined by the majority channel. In this work, we investi- 
gate the contact conductance of single-atom magnetic junc- 
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tions with different chemical species and structures. Our cal- 
culations show that for the Ni-Co/Co(l 1 1) junction the con- 
tact conductance is also mainly determined by the majority 
channel: The majority and minority conductances are 0.38 
and 0.13 Go. respectively. However, for the other three junc- 
tions, the minority conductance is larger than the majority one. 
This indicates that the major contribution to the total conduc- 
tance, from the majority or the minority, is also dependent on 
the junction species and structure. 

In conclusion, by performing first-principles quantum 
transport calculations we have investigated the spin-dependent 
conductance of the four single-atom magnetic junctions. For 
the Ni-Co/Co(l 1 1) junction the conductance is very close 
to Go/2, while for the Co-Co/Co(lll), Ni-Ni/Ni(lll) and 
Ni-Ni/Ni(001) junctions the conductances are far away from 
Go/2 and differ significantly from each other. It was found 
that the spin-up conductance is Uttle influenced by the junc- 



tion species, and therefore, this difference is mainly due to the 
variation of the spin-down component determined by the mi- 
nority dz orbitals which are sensitive to the junction species. 
On the other hand, the junction structure can affect both the 
spin-up and spin-down conductances but not as remarkably as 
does the junction species for the spin-down one. Our calcula- 
tion shows that the previously found Go/2 conductance does 
not generally exist for single-atom magnetic junctions since 
the conductance is largely dependent on the junction species 
and is also influenced considerably by the junction structure. 
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